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ABSTRACT 

We present extragalactic number counts and a lower limit estimate for the 
cosmic infrared background (CIRB) at 15yum from AKARI ultra deep mapping 
of the gravitational lensing cluster Abell 2218. These data are the deepest taken 
by any facility at this wavelength, and uniquely sample the normal galaxy pop- 
ulation. We have de-blended our sources, to resolve photometric confusion, and 
de-lensed our photometry to probe beyond AKARFs blank-field sensitivity. We 
estimate a de-blended 5 a sensitivity of 28.7 fi J j. The resulting 15 /im galaxy 
number counts are a factor of three fainter than previous results, extending to a 
depth of ~ 0.01 mJy and providing a stronger lower limit constraint on the CIRB 
at 15 /im of 1.9±0.5nWm"^sr"^ 

Subject headings: galaxies: clusters: individual (Abell 2218) — galaxies: evolution — 
infrared: galaxies 
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1. Introduction 

The cosmic infrared background (CIRB) is dominated by the dusty-emissions from 
star-forming galaxies, and therefore traces the dust-enshrouded star-formation over the 



history of t he Universe. To inter 
models (e.g. 



Granato et al. 



2004 



jret the CIRB in terms of galaxy formation and evolution 



Pearson 



120051 ) , it is necessary to resolve the monochromatic 
backgrounds into their individua l galaxies. There is a strong correlation between mid-IR 



and far-IR star-forming galaxies fiChary fc Elba2 



2001 



Elbaz et a 



2006 



2002), there 



Devlin et al. 



' ore ga laxies 



20091 ) must 



responsible for the CIRB peak ~140/im-200/im (iDole et al.l 
also dominate the CIRB at shorter wavelengths, i.e., mid-IR <60/im. Recent results from 
Spitzer, SCUBA and BLAST data have shown, via stacking analysis of 24 //m sources at 
longer wavelengths, that 24 /im-selected populations account for the bulk of 70 /im, 160 /im 
and 250 /xm back g rounds and also dominate the 350 Atm, 450 ^tm and 500 /im backgrounds 



(IDole et al. 



2006 



Serjeant et al. 



2008 



Devlin et al. 



20091 ). In contrast, at 850 um. the 



24 um population only resolves around one quarter of the background ( ISerjeant et al 



20081). Previous low er estimates of the CIRB at 15 /im, from lensed ISOCAM data (e.g. 



Metcalfe et al. 



20031 . hereafter MET03), have successfully resolved the contribution from 
galaxies highly luminous in the IR, however the greater depth achieved by these data gives 
a more representative sample of galaxy populations. 



A major challenge for deep IR observations is fl uctuations from confusion noise, which 



presents a fundamental limit to blank-field surveys (jCondon 



19741 ). Exploiting strong 



gravitationa. 



(jSmail et al. 



lensi ng offers a way to probe beyond the inherent blank-field confusion limit 



19971 ). Where lensing increases the apparent surface area of a background field, 
the sources within that area are viewed at a lower number density in comparison to the 
un-lensed situation, and the preservation of surface brightness leads to amplified observed 
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flux densities. Lensing therefore offers a twofold confusion-beating effect. Reduction in 
the observed area leads to source number-density depletion, and both area correction an d 



19951). 



flux correction are required to recover the true galaxy-number counts ( iBroadhurstl 

In addition to faint source confusion below the detection limit, a significant proportion of 

extractions from a conf used image may be blen ds of two or more sources, so photometric 



de-blending is required ( jRodighiero et al. 



20061 ). 



In this Letter we present new 15 /im galaxy number counts and 15 //m integrated 
light (IGL15) estimate. In Section 2 we summarise the AKARI data and data reduction. A 
data analysis description is given in Section 3, and the results are presented in Section 4 
and discussed in Section 5. 

Throughout this Letter we assume fiat ACDM cosmology with f2jv/ = 0.3 and 
i7o = 70kms-^Mpc-^ 

2. Data acquisition and analysis 
2.1. Data 

AKA RI 15 um observations of A22 18 were taken with the L15 filter of AKAUrs IRC 



{AKARI: 



Murakami et al. 



f l2007h : IRC: 



Onaka et al 



f l2007h l The IRC has a wider field 



of view of lO'xlO' in comparison to the Infrared Array Camera (IRAC) aboard Spitzer, 
offering fuller coverage from 2/im-24yum. 19 pointings were acquired with the astronomical 
observation template IRC05, which is designed for deep observations and performs no 
dithering, however a nominal positional offset was applied in-between pointings. 



-5- 
2.2. Data reduction 



The data were reduced using the standard IRC-pipeline, version 20070912 (JLorente et al. 



2008 



hereafter IRC-DUM). During the pipehne each pointing was divided into its constituent 
long and short exposures. The short exposures do not significantly add to the depth or 
quahty of the final frame, so they are discarded. Each long exposure is 16.5 seconds, and 
up to 30 were average-combined per pointing, giving 19 reduced frames. We used optical 
data to register the frames astrometry. 

The frames can be grouped based on the relative scanning direction of the IRC with 
A2218. An interval of roughly six months between the 10*'^ and ll*'^ pointings gave ~180° 
difference in the orientation of the first 10 frames (hereafter L15-A) and the final 9 frames 
(hereafter L15-B). This time interval led to an increase of bad pixels in the L15-B data, 
due to detector degradation. The L15-B data also suffer more severally from scattered 
light, a problem noted in the IRC-DUM. The scattered light is partially addressed by 
the pipeline, but remains an issue for several of the frames. A low-frequency sky noise 
is experienced by all post-pipeline frames. Combining the post-pipeline frames gives a 
significantly uneven background structure, which is detrimental to subsequent photometry. 
We therefore subtracted a median-filtered sky model, generated per frame, using a kernel 
width of 21.5". A comparison of photometry taken for the image combined post-pipeline 
and the image combined after the additional sky-subtraction, showed good agreement at 
the bright end and a systematic shift at the faint end, attributable to the sky-structure 
present in the non-filtered image. We, therefore, concluded that the median-sky subtraction 
removes systematics associated with the extended sky-structure, without detriment to 
source photometry. 

The final frames were average combined, giving an image (hereafter L15-image) 
with total integration time of 8460 seconds, full width at half maximum (FWHM) of the 
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point-spread function (PSF) estimated at 5.96" and a pixel scale of 2.39". Figure 1 shows 
the post-pipeline combined image compared to the L15-image. 



3. Analysis 
3.1. Source extraction 



A 5cr source extraction was performed on the L15-image, with DAOFIND (iStetson 



19871 ). Combining the L15-A frames and the L15-B frames into two 'half-images' gave 
the means for a robust reliability check. Each source was examined individually in the 
'half- images', and those appearing at corresponding coordinates within both images were 
assumed to be real. 



3.2. Sensitivity 

To estimate the map sensitivity we took aperture photometry at random positions, 
excluding the edges. An aperture radius of 5.96" was used, and full flux densities 
were obtained using an aperture correction of 1.44 (See section 3.7) and the IRC-DUM 
ADU-to-/iJy conversion factor of 1.69 /iJy ADU^^. Fitting a Gaussian, with standard 
deviation 8.3 //Jy, to the resulting distribution gave a 5cr sensitivity estimate of 41.7 //Jy. 
For the PSF-fitted catalogue (see Section 3.5) the sensitivity was estimated by comparing 
the input and output photometry for artificial sources introduced to the L15-image, giving 
a 5cr sensitivity of 28.7 /xJy. 
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3.3. Multi-waveband counterparts 



We have multi-waveband coverage of A2218, taken by several facilities: HST WFPC2 
F450, F606 and F814, Palomar 200 in ch Hale u', V, B, i' and WHT's INGRID Ks and J 



(Small et al 



2010 



2001a : 



Ziegler et al 



in preparation); AKARI ^11 ( 



200 Ih: SmtzerlR AC Ch 1 to 4 and MIPS 24 /xm flEgami 



Ko et al. 



20091 ) ■ Figure 2 illustrates the A2218 coverage 



provided by this data set. 

To help identify blended sources within the extracted catalogue, a multi-waveband 
counterpart identification was performed. Potential counterparts for each L15-source 
were identified via a centroid search within a radius of 3.8" (0.64 x FWHM) from the 
L15-centroid. For each source, a comparison of postage-stamp images across the available 
wavebands was performed to identify the main counterpart (brightest) plus subsidiary 
counterparts (less bright) and extra sources in the field not initially identified via the 5 a 
extraction and within a radius of 18.0". 

3.4. Field distortion 



A positional discrepancy between the L15-image and the counterpart images was 
identified during the counterparting process giving a spatially-varying PSF in the L15-image. 
Cubic polynomial coefficients were derived to map the L15-B frames onto the L15-A frame, 
and the resulting frames onto the Palomar and IRAC images. The L15-image was 



recombined and the empirica 



DAOPHOT package (IStetson 



PSF was reconstructed, using PSTSELECT and PSF of the 



19871 ). and showed no spatial variability. The source catalogue 



was then re-centred and the counterparting rerun. 



3.5. De-blending and PSF-fitting 

Simultaneous PSF-fitting was performed on the L15-image using the full post- 
counterparting catalogue, which offers the benefit of positional priors. Constructing a 
reliably representative empirical PSF from a confused image is challenging, so we used 
all suitable sources available to statistically reduce noise in the PSF's tail. The PSF 



was r efined following the iterative method outlined in the DA0PH0T2 manual, (jStetson 



2OOOI ). The PSF radius was set to 17.9", which collects approximately 100% of the flux for 
non-extended sources according to the IRC-DUM and a plot of normalised pixel value as a 
function of radius, for s ources in the L15-image. The res ulting empiric al PSF was used to 



CLEAN (JHogbom 



19871), giving a PSF-fitted 



1974f ) the L15-image with ALLSTAR f lStetsonl 
source catalogue of 918 sources. The increase in the number of sources corresponds to a 
~40% improvement in the completeness (see Section 3.6). 

3.6. Completeness 

Two separate Monte Carlo completeness tests were run to represent the 5a catalogue 
and PSF-fitted catalogue. The first test used the established method of adding randomly- 
placed artificial sources to the L15-image, separated from the 5cr catalogue, then performing 
an extraction on the results followed by aperture photometry of the extracted sources. The 
artificial sources were randomly scaled within defined flux density bins covering the source 
catalogue's flux range. For the first test the artificial source positions were generated with 
sufficient separation to avoid self- confusion. The 25" minimum separation was derived by 
plotting normahsed pixel values as a function of radius for bright well-separated sources. 
Twice the radius where the median pixel values disappear into the background was chosen. 
This test was repeated until around 20,000 sources per bin were achieved. The first test 
was adapted to represent the PSF-fitting of a photometrically-confused environment. Input 
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positions were still generated randomly and kept at a distance from known source positions, 
although this limiting separation was reduced to 19". A self-separation was imposed, 
but only to reject equal random positions. To reflect the use of positional priors and the 
re-centring carried out by ALLSTAR, the randomly-generated input positions were used as 
the ALLSTAR input rather than the extracted positions. This second test was repeated 
until around 30,000 sources per bin were achieved. 

Completeness was defined as the fraction of recovered sources per bin. The results 
of the first tests show that the L15-image is 10%, 50% and 90% complete down to 20.2 /iJy, 
30.7 /iJy and 46.8 /iJy, respectively. For the second test, the L15-image is 10%, 50% and 
90% complete down to 12.2 /iJy, 20.0 /iJy and 31.5 /iJy, respectively. 



3.7. Multi- waveband photometry 



Small et al. 



(12001bl ). 



HST photometry was obtained from the published catalogue of 
IRAC aperture photometry was taken with an aperture radius of 2.44" and an annulus of 
radii 14.6" and 24.4", and the published IRAC aperture corrections were applied. For the 
remaining counterpart images, ape rture photometry was taken and a growth-curve aperture 



correction method was employed (JHowell 



1989 



Stetson 



1990l ). For each image a median 



growth-curve was empirically constructed using aperture photometry taken for bright and 
well-separated sources, with concentric apertures of increasing radii. The u' to Ks images 
were better represented by two growth-curves, one for point-like sources and the other for 
elliptical sources, which are not significantly extended. Aperture corrections were chosen 
on a source-by-source basis to minimise contamination from neighbours. For the L15-image 
an aperture correction for a radius of 5.96" was derived, using the empirical PSF and a 
comparison of the PSF-fitted photometry and aperture photometry. 
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3.8. Photometric redshifts 

Two codes that utilise a minimum x^ spectral energy distribution (SED) fitting method 
were applied to estimate photometric redshifts for L 15-sources with photo metry coverage 



20081 ) is suitable for 



in four or more filters, shortwards of 11 /im. EaZy (iBrammer et al.l 
data sets with few or biased spectroscopic redshifts (z^pec), such as the z^pec available for the 
LI 5- catalogue, which are mainly biased at the cluster redshift of 0.18. The EaZy theoretical 
SED templates are based on semi-analytical models, and a linear combination of templates 
can be fitted simultaneously. IRAC photometry was included due to EaZy's ability to fit 
photometry up to IRAC CH4, however this is dependent on redshift. EaZy gives the option 
to apply priors aimed at breaking the template colour degeneracies seen with increasing 



Negrello et all (12009 



redshift. Our spectra were also fitted using the photometric code of 
hereafter N09). This code is uniquely optimised for fitting mid-to- far- infrared polycylic 
aromatic hydrocarbon (PAH) and silicate features seen in starburst SEDs. For sources with 
strong mid-infrared PAH features, the comparison of N09 and EaZy redshifts was consistent 
with a slope of 1. A robust catalogue of photometric redshifts was constructed using a visual 
triple-check per source to reject unreliable estimates. The best SED fits from EaZy and 
N09 and the source morphology were visually compared, in context of the redshift estimate 
and probability of the minimum x^. The redshift catalogue was constructed primarily from 
EaZy estimates. For sources with pronounced mid-infrared features and reliable EaZy 
and N09 estimates, not in agreement within their la errors, the N09 estimates were used 
when clearly providing additional constraint from fitting to mid-IR photometry. Cluster 
members were identified from spectroscopic redshifts or during the triple comparison, 
from their typical SED and elliptical morphology. Cluster members represent 16% of the 
total catalogue, including all significantly extended sources. All cluster members were 
subsequently removed. 31% of the remaining catalogue are without a redshift estimate, 
either due to a lack of multi-wavelength coverage or unreliable photometric estimate. For 
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these sources a redshift of 1.04 ±0.67 was assigned, which is the median of the redshift 
catalogue with la errors. Substituting a value of 2.0 or 3.0, in place of the median redshift 
value, gave no significant difference for the resulting number counts. 

4. Results 



4.1. A2218 mass model 



20071 ) 



Magnification corrections (/i) were obtained using LENSTOOL (iJuUo et al. 
which required as input a mass model of A2218 and the positions and redshifts for all 
sources beyond the cluster distance. A pseudo isothermal elliptical mass model is assumed 
for the mass distribution of the A2218 cluster members. Strongly lensed arcs and arclets 



are used to constrain the mass distribution and total mass of the cluster ( Eliasdottir et al. 



2007 



Knei 



3 et al 



19961). Spe ctroscopy of arclets has been used to test this model for 



reliability (JEbbels et al. 



19981 ). 



4.2. Galaxy number counts 



Flux densities (S) were corrected prior to counting as S^ 



true 



■-^obs 



Corrections for 



depletion and incompleteness were applied to individual sources during counting, assuming 



the relation ntr 



C(S'obs) 



, where ntme is the true number of sources and the completeness 



((7) is a function of S'obs (rather than 5'true)- De- lensed number counts over bin dS are 
then obtained as ^ = Sntruc- The amplification (-) distribution ranges from 1.0 to 24 
and has a median of 1.2, which reflects the wide area of A2218 covered and the decrease 
of amplification as a function of radius from the centre of the core. The fi distribution 
obtained was not found to change significantly with the variation of redshifts within the 
L15-redshift distribution. 
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Figure 3 shows our Euclidean-normalised differential number counts, in comparison 
to a compilation of previous work and predictions based on galaxy evolution models. The 
median completeness per bin is 100% down to the faintest three bins, which have median 
completeness corrections of 26%, 70% and 96% respectively. Our L15-counts extend the 
faint end of obser ved counts down to ^ 0.01 mJy, which is a factor of 3 fainter in comparison 



to the ISOCAM fICesarskv et al. 



19961 ) lensing survey counts of MET03. Below 0.2 mJy the 



L1 5-counts present a steep sub-Euclidean slope of -1.6, which agrees with the faint slope 



of 



Elbaz et al. 



( 119991 ). The no-evolution model is strongly excluded by all available data 
and there is a general consensus on a significant evolutionary bump, which peaks around 
0.2-0.4 mJy. The comparably steep slope of the L15-counts brighter than the 'bump' is 
the result of de-blending. The Pearson 2007/2010 model (hereafter PIO) predicts that the 
populations dominating the 15 /im counts 'bump' are starbursts (L < 10^^ Lq) and luminous 
infrared galaxies (LIRG; 10^^ L0 <L<1O^^L0), with redshift distributions peaking at 
z = 0.5 and z = 1.2, respectively. At a mean redshift of 0.8, in the bump, S15 = 0.3 mJy 
corresponds to ~4xlO^^ Lq for an M82 SED. 

4.3. Bootstrapping 

Confidence intervals for the differential counts were derived by bootstrapping within 
the photometric and redshift errors, as these are the dominant source of uncertainties for 
the counts. The L15-flux population was re-sampled without bias for the lensed-source 
catalogue. Each sample was randomly assigned flux densities and redshifts within the 
respective 3cr errors, and de-lensed with recalculated magnification corrections. 30,000 
re-sampled populations were generated and differential counts were taken for each. 
Confidence intervals were calculated using the median and standard deviation of the 
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resampled counts. The grey shaded region of Figure 3 shows the resulting 95% confidence 
interval. The divergence of the bootstrapping below 0.02 mJy shows that the upturn of the 
counts in the final bin is not statistically significant. The differential number counts and 
the bootstrapped standard deviation are presented in Table 1. 

4.4. 15 fim integrated galaxy light 



The IGLi5, a lower limit for the CIRB15, can be obtained by integrating the fiux 
per unit area for the corresponding monochromatic number counts. The differential 



contribution to the IGL15 is given by 
dlGL dN , S , 



dS 



dS' 10^0' 



(e.g. 



Elbazetal 



20021) 



where 1 mJy = -^ nW m ^ sr ^ and z^is is the frequency of the 15 yum photons. 



Usi ng the counts of 



mm), 



Gruppioni et al. 



Serjeant et al 



mm; 



Rush et al. 



(I2OOOI re-c alibrated following 



Vaisanen et al. 



(119931 ) and the L15-counts, giving a flux 
range of 0.01 mJy-10,000 mJy, we estimate IGL15 = 1.9 ± 0.5 nWm^^ sr^^. The lensed flelds 
observed by ISOCAM (including A2218) produced 15 //m counts down to 0.03 mJy and 
a lower-limit estimation for the IGL15 of 2.7 ± 0.62 nW m^^ sr^^ (MET03, and references 
therein). These estimates agree, within the errors. Using the lower flux limit of MET03, and 
our methodology, gives an IGL15 of 1.6 ± 0.38 nW m~^ sr~^, which is marginally consistent 
(~2cr) with the MET03 estimate. The PIO model provides an excellent flt to the whole of 
our data. Integrating over the full flux range of the PIO model gives a predicted IGL15 of 
2.3nWm^^ sr^^. If we assume the shape (but not the normalisation) of the PIO counts, 
we derive a slightly better estimate of IGL15 = 2.0±0.4nWm^^ sr^^ at >0.01mJy (see 
Figure 4). In comparison to the MET03 result, which resolved ~ 70% of the CIRB15 into 
individual galaxies, we are x3 deeper and resolve 87% ± 13% and, whereas the Infrared 
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Space Observatory (ISO) surveys mainly sample galaxies with luminosities > LIRG, we are 
probing the more normal galaxy populations. 

5. Discussion and conclusions 

From our de-blended and de-lensed 15 yum counts we have derived an IGL15 estimate 
of 2.0 ±0.4nWm~^ sr~^, down to ~0.01mJy. We conclude that, with respect to the 
PIO model, the AKARI 15 /im data are consistent with having resolved the whole of the 
predicted IGL15. Assuming no radical change between the IR SEP of high-re dshift galaxies 



and those resolved at 15 /im with median redshift of 1.0 (JElbaz et al. 



20021 ) ■ then the 



galaxies resolved by these data represent the bulk of galaxies dominating CIRB peak. 

Figure 4 suggests that in order to resolve 100% of the CIRB15, future observations need 
to probe depths in the region of 1 magnitude fainter than the sensitivity limit achieved by 
this survey, down to at least 5*15 = 1 /iJy. The fir st possible direct measurement constr aints 



of the CIRB15 will come from JWST or SPICA (.Gardner et al. 



2006 



Nakagawa 



200J). 



15 /im stacking analysis of Herschel/ SPIRE and PACS A2218 data, will address how 
representative 15 /im selected galaxies are of the galaxy populations responsible for the 
CIRB at its peak. 
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Fig. 1. — Comparison between the L15- image combined post-pipeline (left) and the L15- 
image combined after the further reduction (right). The same pixel scaling was used to plot 
both images. The low-frequency background structure evident in the post-pipeline image is 
successfully removed by the additional reduction. 
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Fig. 2. — L15-image overlaid with contours indicating the multi- wavelength coverage of 
A2218. 
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Table 1. Lensed ( .L "'^ ) and De-lensed ( .L *''"° ) Differential Number Counts, corrected 

for Incompleteness, and the Associated Bootstrapped Median Number Counts and 

Standard Deviation for the De-lensed Counts. 



Jlll^ l^iiiojy 


dSob >-"'•'> ^-B ; 


<i5t„„ V"'-'> •^-B 


l.OOE-02 


(4.37±0.98)E+06 


(3.06±0.99)E+06 


1.59E-02 


(8.60±0.58)E+05 


(7.44±0.92)E+05 


2.51E-02 


(3.61±0.81)E+05 


(4.63±0.35)E+05 


3.98E-02 


(2.23±0.12)E+05 


(2.02±0.13)E+05 


6.31E-02 


(9.42±9.47)E+04 


(1.08±0.14)E+05 


l.OOE-01 


(4.44±0.52)E+04 


(5.62±0.81)E+04 


1.59E-01 


(1.93±0.27)E+04 


(1.90±0.38)E+04 


2.51E-01 


(7.17±1.31)E+03 


(6.81±1.79)E+03 


3.98E-01 


(2.56±0.62)E+03 


(1.62±0.69)E+03 


6.31E-01 


(2.47±1.23)E+02 


(1.37±1.29)E+02 


1.20E+00 


(0.00±0.00)E+00 


(4.68±6.35)E+01 



^i'truc bootstrap 



(mjy ^dog 2) Crbootstrap(m-Iy ^dog ^) 



2.76E+06 
7.71E+05 
4.27E+05 
2.16E+05 
1.06E+05 
4.73E+04 
1.73E+04 
5.92E+03 
2.24E+03 
2.65E+02 
4.39E+01 



3.54E+05 
7.13E+04 
3.28E+04 
1.80E+04 
1.03E+04 
5.70E+03 
2.73E+03 
1.22E+03 
6.00E+02 
1.37E+02 
4.01E+01 



Note. — Lower bin limits are given 



